Fragile granular jamming 
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Abstract 

^ , We demonstrate experimentally that the route to a jammed state for a set of bi-dispersed frictional disks, 

^ I subjected to uni-axial compression from a random initial unjammed state, consists of a consolidation state, 

■ a fragile jammed state, and finally a rigid jammed state. In the consolidation regime, the pressure on the 

I . sides increases very slowly with the packing fraction (p, and there are no detectable stress chains. In the 

C I fragile jammed state, stress chains are visible, the pressure increases exponentially with (f), and the fraction 

^ ' of moving disks drops exponentially. Eventually, a final regime where particle displacements are below our 

resolution and the pressure varies approximately linearly with is reached. We argue that this scenario is 
OO , generic for athermal frictional compressed particles. 
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When loose granular matter is progressively compressed to higher packing fractions 0, the 
material stiffens at a critical packing fraction cpc owing to development of inter-particle constraints 
- the system is then said to be jammed. This transition from an athermal, particulate gas to a 
disordered solid has been related to the behavior of a wide variety of materials including glasses, 
foams, and colloids in a unifying phase diagram proposed by Liu and Nagel yj, |20. From the 
perspective of real granular materials, however, the effects of inter-particle friction on the jamming 
scenario have not been widely considered yfl. On the contrary, to allow better comparison with 
the frictionless model, various empirical protocols such as mechanical shaking, repeated cycling, 
or high-pressure compaction (useful for relatively soft deformable particles) have been introduced 
in experimental systems to help break down weak frictional contacts. Nevertheless, friction can 
play a significant role in such materials, many of which are technologically relevant (as in the 
compression and sintering of ceramic [Q]). It is thus important to understand the full impact of 
frictional contacts without resorting to mechanical vibration or repeated cycling whose effects are 
difficult to quantify. 

Experiments [|5i|6D have confirmed many features of the jamming transition, such as power-law 
scaling of the pressure P and of the mean coordination number Z with increasing [[SO, and a 
maximum in the pair-correlation function [1^. Given the unquantified mechanical vibration, cy- 
cling or highly-plastic deformation employed in the experimental work, an analysis of results that 
are strictly athermal, are weakly deformed during compression, and that fully manifest frictional 
effects is called for. 

In this Letter, we report a quasi- static experiment across the jamming threshold, performed on 
a two-dimensional bi-disperse system of photo-elastic disks. We show that the pressure scaling 
is remarkably different between the first and subsequent jamming cycles. In the first cycle as 
the system's boundaries are moved in to achieve an increasingly tighter packing, the pressure 
exhibits a short ramp followed by an exponential rise P oc e"^/^^. Eventually the pressure deviates 
from exponential behavior and scales linearly with respect to the packing fraction P oc 0. When 
decompressed and compressed again, however, the intermediate region with exponential scaling 
disappears and one recovers the approximate linear scaling [fiol] . The exponential increase in 
the pressure is reflected in an exponential decrease of the mean and variance of the domain size 
enclosed by stress chains above a certain cutoff and in an exponential decrease of the variance of 
particle displacements and of the fraction of particles moving more than about 1% of a particle 
diameter. Our measurements demonstrate that there are three regimes of frictional jamming during 
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FIG. 1: Schematic of the experimental setup. The system consists of 950 large and 950 small disks (ratio of 
radii 1:1.5). Two movable boundaries at opposite ends control the system's packing fraction. Force sensors 
labeled A through F in the schematic measure the boundary pressure. The image is contrast enhanced data 
for (p = 0.8113. 

the first lateral compression: a consolidation regime where particles move freely; a fragile jammed 
region where pressure increases exponentially, particle displacements become smaller and less 
frequent, and a stress chain network forms; and finally a fully jammed state with approximately 
linear dependence of pressure on area fraction, very small displacements of bulk clusters, and an 
approximately static stress chain network. Measurements were also performed on disks machined 
with varying friction coefficients (Lexan, Teflon, and photo-elastic disks coated with graphite 
dust). The intermediate regime with an exponential rise in P was observed in all measurements 
over the same range of pressure, with exponents xp varying with the material bulk modulus. These 
results will be discussed in detail elsewhere. 

Experimental set-up — Figure [T] shows a schematic of the experimental setup. The system 
consists of a bi-disperse mixture of 950 large (diameter D = 0.9525 ± 0.0025 cm) and 950 small 
(diameter d = 0.635 ± 0.0025 cm) polymer disks of thickness 0.508 cm that are bi-refringent 
under stress (photo-elastic) and that have an inter-particle friction coefficient of about 0.2. The 
disks are placed in a chamber, with dimensions L = 25.9 cm and W = 48.2 cm where L is the 
compression direction, consisting of two glass plates held 0.635 cm apart by means of an acrylic 
frame that runs along the system's perimeter. Two movable boundaries are placed on the acrylic 
frame with aluminum plates that can slide back-and-forth within the chamber from opposite ends. 
The transverse boundaries are held fixed. The positions of the movable boundaries is controlled by 
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two micrometers with a precision of 0.001 cm. Taking variations in radii into account for the given 
system of 1900 disks, this translates to a precision A0 = 1 x 10~^ in the packing fraction, and 
hence serves as the lower bound on the quasi-static step-size (A0). All measurements reported 
in this letter, however, were made at a quasi-static step-size of A0 = 1 x 10^"^. The packing 
fraction (0) is defined as the ratio of the area occupied by the disks to the total chamber area. 
The packing fraction is therefore controlled by changing the chamber area in this experiment. A 
set of six sensors (labeled A through F in Fig.[T]) placed along the boundaries measure the global 
pressure. The noise floor of these sensors is of the order of 0.1 N/m. Along the compression 
direction, the sensors feel the friction of the movable boundaries with the glass bottom, which 
was measured to be about 1 1 N/m, accordingly this serves as the true zero of pressure measured 
along the compression axis. Visual measurements using a Nikon D-90 camera (12.3 megapixel 
resolution) yield positions and displacements of individual disks as well as the stresses inferred 
from the photo-elastic response of the disks ||5|]. We determined the local disk contacts having a 
stress above our experimental light intensity threshold (corresponding to about 1 N) and satisfying 
local force balance constraints. The entire stress chain network satisfies physical constraints and 
provides a good characterization of the distribution of local forces to complement the pressure 
measurements. 

Protocol — We initially place all disks in the interrogation chamber in random positions and 
move the boundary plates to an initial packing fraction chosen a priori to fall comfortably in 
an unjammed state. The disks are subject to friction (both with the bottom boundary and with 
each other) and the system jams below the random close packed density (prep ~ 0.84. After each 
quasi-static step, a 10 s time-trace of all six boundary sensors is collected at a sampling frequency 
of 1 KHz followed by a digital image of the whole system. The time-averaged value of the trace 
constitutes the pressure measured by the boundary sensor at a particular value of 0. The boundaries 
are then moved through a quasi-static step (A0 = 1 x 10"''), and the procedure is repeated. We 
do not tap/anneal the system after each quasi-static step — as is often done in numerics — nor do 
we excite the system vibrationally to breakdown frictional effects: we are interested in studying 
the system evolution at "zero temperature" as it occurs under the control of a single parameter, 
namely (p. 

Three regimes — The argument for a new, intermediate, fragile state is best appreciated by 
comparing the pressure signal for the first and second jamming cycles. Figure [Ja) plots P from 
one compression boundary sensor versus (p for the first and second jamming cycles. The con- 
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FIG. 2: (Color online) (a) P vs (p - linear scale - for the first (solid circles) and second (solid squares) 
jamming cycles. The pressure scaling is gradual for the first cycle as compared to the more abrupt transition 
during the second cycle where (j)c is indicated. The inset plots the pressure from all four boundaries (same 
vertical and horizontal scale). No anisotropy is observed in the pressure signal, (b) P vs for the first 
jamming cycle - log-lineai^ scale - shows the existence of an intermediate regime where pressure scales 
exponentially. The solid line is an exponential fit to the data P oc e^^^'' with xp = 0.002, and the dashed 
line is a linear fit. Values of 0i and 02 are indicated in the plot. 

tinuous increase in P for the first cycle is qualitatively different from the more abrupt change in 
slope for the second cycle. The lateral shift in P is a signature of a friction-induced hysteretic 
response which progressively shifts (pc to higher values as the system is repeatedly jammed and 
unjammed (discussed in detail elsewhere). The vertical shift in P in the unjammed regime for 
the first cycle is traced to the friction between the movable boundary plates and the glass bottom. 
During the first cycle, the boundaries move in and push the disks towards a jammed configuration. 
When unjammed and jammed again, the contacts that were developed at the end of the first jam- 
ming cycle are re-activated, and the stresses build up as the system is subjected again to uni-axial 



5 




FIG. 3: (Color Online) Mean (•)and variance (o) of stress chain domain area normalized by total area, am, 
(left axis) or by the mean area formed by connections between different combinations of 3 disks in a close 
packed triangular array, am, (right axis), respectively. Insets show stress networks and corresponding 
domains for (p = 0.8104 (left) and (f) = 0.8148 (right). The dashed lines are fits to a combined hnear 
dependence with a decaying exponential (see main text). 

compression. 

In order to better understand the smooth increase in P for the first jamming cycle, Fig. [2b) 
plots P versus on a log-linear scale. One sees three distinct regimes in the pressure curve. In 
the unjammed (consolidation regime) P is essentially flat, modulo a shallow ramp due to friction 
between disks and bottom glass plates. The second regime is characterized by an exponential 
increase in P beginning at (pi ^ 0.8093. The solid line in Fig.[2tb) is an exponential fit to the 
data P oc e*^/^^ with xp = 0.002. The pressure eventually deviates from this exponential regime 
and settles to a linear scaling P oc at 02 = 0.8124. The inset in Fig. [2] (a) plots P from all four 
boundaries against on a log-linear scale. The remarkably good collapse of the pressure curves 
atop one another strongly suggests that anisotropic effects arising from uni-axial compression are 
not detected by the pressure sensors. This isotropy may result from our geometry for which the 
aspect ratio is L/ 14^ = 0.54 and would probably not persist for large aspect ratios, i.e., L/H^ >> 1. 

This exponential pressure scaling for the first jamming cycle contradicts the predicted O, lZ|, [SQ 
pow er-law for P across the jamming transition. This power-law scaling (approximately linear 
lIlOll ) is, however, recovered for subsequent jamming cycles as shown for the second jamming 
cycle in Fig. [2k where (pc = 0.8113 is determined by the point at which P starts to rise. 

To elucidate the nature of the different regimes that we observe in global pressure measure- 
ments, we consider local properties of disk configurations, namely the structure of stress chains 
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and the measurements of disk displacements. The identified stress chains enclose domains with 
no measured stress as illustrated in the insets of Fig.[3l We characterize the stress chain networks 
by the mean and variance of the fractional domain area (relative to the total cell area on the 
left axis) and also by the mean and variance of the domain area but now normalized (right 
axis) by the average of the area enclosed by a triangular arrangement of small and large disks (all 
large, all small, two large, two small). As indicated in the inset images, the mean and variance 
decrease rapidly between 0.810 (the lowest value of cf) for which we could determine the stress 
chain network) and 02 = 0.8124. For higher > 0.8124, the mean and variance decrease linearly 
with a small slope. Fits to an exponentially decreasing function e~'^^^^ as shown in Fig. [3] yields 
the same value of Xs = 0.00035. The exponential decrease in and is another signature of 
the fragile jammed state. 

We now consider the behavior of disk displacements. Figure [4^1 shows the fraction of disks 
moving a distance greater than about 1 % of a mean disk diameter Iil2il as a function of where 
is the number moving and Nt = 1900 is the total number of disks. The fraction is constant at about 
0.3 up to 01 = 0.8094, after which it decreases rapidly up to the jamming value 02 = 0.8123. The 
decrease is consistent with an exponential e~'^^^'^ with xn ~ 0.001. The variance of individual 
disk displacements a (normalized by = (D + d)/2 = 0.794 cm) relative to a state near the 
jamming threshold lllin is shown in Fig. |4] as a function of 0. One again observes two distinct 
regimes: a linear one that corresponds to the unjammed (consolidation) regime where the pressure 
curve in Figure [2tb) is flat, and a second one in which the variance decreases exponentially with 
a oc e"*/^"^ with Xa = 0.0007. This exponential drop in displacement variance occurs over the 
same interval in as the exponential regime of the pressure curve in Fig. [2tb) as indicated by 
the arrows depicting 0i and 02. The presence of disk displacements in the exponential regime, 
where a percolating force network already exists, suggests that the particles that are part of the 
force network still undergo small displacements and deformations, which in turn allows visible 
displacement inside the network region, eventually leading to the refining of the network. This 
picture is reinforced by the data in the inset of Fig. Il^b) where the difference in the force chain 
network between = 0.8095 and = 0.8105 is shown as black lines and the spatial distribution 
of the magnitude of particle displacements is also shown: The correlation of new stress chain 
creation (the differences) and the particle displacements is striking. 

Figures |2] and IH when considered together, present a picture of a fragile system where part of 
the system is jammed (as evidenced by an exponential increase in pressure) whereas another part 
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FIG. 4: (Color online) (a) Fraction of disks moving Nm/NT vs (p showing constant region for < 0i and 
exponential decay for (j)i < (f) < with xn ~ 0.001. (b) Normalized displacement variance cr/ti relative 
to a state near 02 showing a linear decrease for (j) < and an exponential decrease with Xa = 0.0007 for 
(t)i < (j) < 4>2- (c) Image of the superposition of the difference in the stress network (black lines) and the 
magnitude of disk displacements between cj) = 0.8095 and (p = 0.8105 Aitows indicate 0i and (j)2- 



is still floppy (as seen from the variance of disk displacements, albeit decreasing exponentially 
with (j), and the corresponding decrease in the mean and variance of domain area in Fig.[3]). These 
results collectively form the crux of the argument presented in this Letter arguing in favor of a 
fragile jammed state that separates the unjammed and the rigid one. 

Discussions — We have presented compelling evidence for the existence of an intermediate 
fragile state separating the unjammed consolidation state from the rigid jammed state of athermal 
particles with frictional contacts. The fragile state is characterized by an exponential rise in pres- 
sure against a power-law scaling predicted by field theories [|8|], simulations [|7D and experiments 
that minimize friction through external perturbations or cycling [5]. Friction seems to play an 
essential role in providing the mechanical support leading to a marginally stable, yet fragile solid. 
These findings may be related to the experiments of Ref.[161] where the pressure shows a similar 
smooth rise in a very different system, suggesting that our scenario is a generic one for athermal 
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frictional compressed particles. 

Comparing the exponential behaviors of different quantities in the fragile jammed state pro- 
vides an interesting but incomplete part of our work. Since Xs = 0.00035 measures the behavior 
of domain area, a domain length ^ = ^/a^ would be expected to have twice the decay constant 
(owing to the square root) or 0.0007. This is just the value of Xa- One might argue that the global 
pressure is some average over the network stress chain, say P ~ J^^sc/^^, which implies that 
forces would increase exponentially with xf ~ 0.002 — 0.0007 = 0.0013. It would also be inter- 
esting to see if there is any relation between the intermediate regime and the so-called intermediate 
phase that occurs in rigidity percolation [9]. In both cases, a direct and precise measurement of 
stress chain magnitudes would allow a quantitative test. In any event, the robust exponential be- 
havior observed here experimentally certainly calls out for a theoretical explanation. 

This work was carried out under the auspices of the National Nuclear Security Administration 
of the U.S. Department of Energy at Los Alamos National Laboratory under Contract No. DE- 
AC52-06NA25396. 
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